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Adsorption of P-H2 films on Alkali metals substrates at low temperature is studied theoretically 
by means of Path Integral Monte Carlo simulations. Realistic potentials are utilized to model the 
interaction between two P-H2 molecules, as well as between a P-H2 molecule and the substrate, 
assumed smooth. Results show that adsorption of p-Hb on a Lithium substrate, the most attractive 
among the Alkali, occurs through completion of successive solid adlayers. Each layer has a two- 
dimensional density 9 e ~ 0.070 A~ 2 . A solid p-Hh monolayer displays a higher degree of confinement, 
in the direction perpendicular to the substrate, than a monolayer Helium film, and has a melting 
temperature of about 6.5 K. The other Alkali substrates are not attractive enough to be wetted by 
H2 at low temperature. No evidence of a possible superfluid phase of P-H2 is seen in these systems. 

PACS numbers: 67.70+n, 61.30.Hn, 68.08.Bc 67. 



I. INTRODUCTION 

Much theoretical and experimental work has been de- 
voted to the investigation of films of highly quantal flu- 
ids (e.g., condensed Helium) adsorbed on a variety of 
substrates, from the most strongly attractive, such as 
graphite, to relatively weak ones, such as Alkali metals 1 . 
The motivation underlying this effort is the understand- 
ing of the fascinating properties that quantum films dis- 
play, often quite different than those of bulk materials. 

Thin films of para- Hydrogen (P-H2) are of consider- 
able interest. One of the reasons is that a fluid of p- 
H2 molecules is regarded as a potential "second super- 
fluid", due to the light mass and the bosonic character 
of its constituents 2 . In bulk P-H2, superfluidity is not 
observed, because, unlike Helium, H2 solidifies at a tem- 
perature (T 13 K) significantly higher than that (~ 4 
K) at which phenomena such as Bose Condensation and, 
possibly, superfluidity (SF) might occur. This is due the 
depth of the attractive well of the potential between two 
Hydrogen molecules, significantly greater than that be- 
tween two Helium atoms. Several, attempts have been 
made 3-4 - 5-6 to supercool bulk liquid p-B.2- but the search 
for SF (in the bulk) has so far not met with success. 

Confinement, and reduction of dimensionality, are 
widely regarded as plausible avenues to the stabilization 
of a liquid phase of P-H2 at temperatures sufficiently low 
that a SF transition may be observed. Indeed, com- 
puter simulations yielded evidence of superfluid behav- 
ior in very small (less than 20 molecules) p-H.2 clusters^, 
and claims have been made of its actual experimental 
observation^. Theoretically, SF has also been predicted 
to occur in a strictly two-dimensional (2d) P-H2 fluid em- 
bedded in a crystalline matrix of Potassium atoms 9 . 

Interesting as the above predictions undoubtedly are, 
they pertain to physical systems that are not easily re- 
alized in a laboratory, nor are controlled measurements 
of their properties always practical. On the other hand, 
adsorbed films of P-H2 are readily accessible experimen- 
tally, and in fact have been extensively investigated, on 
different substrates. For example, the phase diagram and 
structure of monolayer P-H2 films adsorbed on graphite 
have been studied by various technique9i2*ii*i2ii£. One 



of the most remarkable aspectsi 3 - is that the melting 
temperature T m of a solid P-H2 monolayer can be sig- 
nificantly less than bulk p-H2- The question arises of 
whether a considerable reduction of T m could be achieved 
on a weak substrate, such as that of a Alkali metal. It 
has been shownii that, for a 4 He monolayer adsorbed on 
Lithium, zero-point motion in the direction perpendicu- 
lar to the substrate is quite significant, accounting for a 
doubling of the kinetic energy per atom, with respect to 
strictly 2d He. Such zero-point motion, conceivably im- 
portant in a P-H2 film as well, could result in an effective 
screening of the interaction of P-H2 molecules, leading to 
a suppression of the melting temperature, perhaps to the 
point where novel quantum many-body phenomena may 
become observable in the liquid phase. 

In this work, attention is directed to the physics of 
P-H2 films adsorbed on Alkali metals substrates. Be- 
sides the issue of SF, this system has elicited interest 
among condensed matter scientists for its intriguing wet- 
ting properties. It was suggested^, and experimentally 
verifiedi 5 ^iii*i£, that H2, like 4 He, ought not wet Cesium 
(Cs) or Rubidium (Rb) substrates at low temperature, 
due to the relatively shallow depth of the substrate- 
adatom potential. It is not known whether a thermo- 
dynamically stable P-H2 film should exist, at low T, on 
any of the other Alkali substrates, namely Potassium (K), 
Sodium (Na) or Lithium (Li). Particularly intriguing is 
the question of the existence of a P-H2 monolayer on 
these substrates, and what its melting temperature T m 
should be. 

Aside from the study of the phase diagram of 2d p- 
H2 by Gordillo and Ceperlej*^, and of a P-H2 surface 
by Wagner and CeperlejiiSii^, theoretical calculations 
for realistic models of adsorbed P-K2 films have focused 
on graphite^ .. Only recently have the first theoreti- 
cal studies of P-H2 films on Alkali substrates been car- 
ried out, based on Path Integral Monte Carlo (PIMC) 
simulations^^ as well as on semi-empirical Density 
Functional methods 2 ^,. PIMC studies have shown evi- 
dence of a wetting transition for liquid P-H2 on Rb and 
Cs at temperatures well above 20 K. 

Here, we present detailed PIMC results for a model 
of P-H2 film adsorbed on Li and Na substrates, in the 
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temperature range between 2 and 13 K, i.e., significantly 
lower than that studied in Ref. 0. The most accurate 
potentials currently available are utilized to describe the 
interactions among P-H2 molecules, as well as between 
the molecules and the substrates, which are regarded as 
smooth (i.e., corrugation is neglected). 

The main results of this study are the following: 

1. No stable p-Eb film forms on a Na substrate, in 
the T — » limit; the same is expected to hold true 
a fortiori on a K substrate, which is more weakly 
attractive than Na. 

2. Adsorption on a Li substrate occurs through com- 
pletion of successive solid layers. A stable solid 
P-H2 monolayer exists, whose equilibrium coverage 
(i.e., 2d density) is 6 e = 0.070 ± 0.01 A~ 2 . This 
is the same value (within statistical uncertainty) 
found on graphite^ ., only slightly greater than the 
theoretically computed equilibrium density of 2d p- 
H 2 (see Ref. i). 

3. The spread of a p-Eb monolayer in the direction 
perpendicular to the substrate is significantly less 
than for a 4 He monolayer, and its melting temper- 
ature is ~ 6.5 K, i.e., very close to what has been 
theoretically predicted^ for 2d P-H2 ■ 

The quantitative similarity of the results obtained in this 
work to those involving different substrates, of greatly 
varying strength and corrugation, suggests that the basic 
physics of this system is driven by the interaction among 
hydrogen molecules, and that the role of the substrate is 
rather marginal. 

Although the PIMC method utilized in this work al- 
lows for the sampling of permutations of particles, which 
is essential in order to reproduce in the simulation any 
effect due to quantum statistics, permutations are not 
seen to occur, in the temperature range explored. This 
is consistent with the conclusion reached in previous work 
by others^. In the crystal phase, permutations are sup- 
pressed (as in most solids) by the localization of p-Eb 
molecules; at higher temperature, on the other hand, 
though the crystal melts and molecules are less local- 
ized, they also behave more classically, as their thermal 
wave length decreases. Consistently with permutations 
not being important, i.e. P-H2 molecules obeying Boltz- 
mann statistics in the temperature range explored here, 
no evidence of SF can be seenS. 

The remainder of this manuscript is organized as fol- 
lows: in the next section, the model Hamiltonian is in- 
troduced, and a brief description is provided of the most 
important aspects of the computational method utilized 
(PIMC), of which comprehensive reviews exist in the lit- 
erature; we then illustrate our results and outline our 
conclusions. 



II. MODEL 

Consistently with other theoretical studies, our sys- 
tem of interest is modeled as an ensemble of N P-H2 



molecules, regarded as point particles, moving in the 
presence of an infinite, smooth planar substrate (posi- 
tioned at z = 0). The quantum- mechanical many-body 
Hamiltonian is the following: 

*2 N N 
i—l i<j i—1 

The system is enclosed in a vessel shaped as a par- 
allelepid, with periodic boundary conditions in all di- 
rections. Here, m is the P-H2 molecular mass, V is 
the potential describing the interaction between two p- 
H 2 molecules, only depending on their relative distance, 
whereas U is the potential describing the interaction of a 
helium atom with the substrate, also depending only on 
the distance of the atom from the substrate. The Silvera- 
Goldman potential 2 ^ was chosen to model the interac- 
tion between two P-H2 molecules; this potential has been 
shown to provide an acceptable quantitative description 
of bulk P-H2 in the condensed phase 25 . 

The following potential, proposed by Chizmeshya, 
Cole, and Zaremba»Si was adopted to describe the inter- 
action of a P-H2 molecule with a smooth substrate (i.e., 
the U term in ) 

jTi \ tt f-t , \ -ax Cvdw f2(P{z)(z — Z v d w )) . , 

U[z) = U {1 + az)e -—— -3 (2) 

\Z Z vc [w )' 

with/ 2 (x) = 1- e~ x (l + x + x 2 /2) and (3(x) = a 2 x/{l + 
ax). The first term in @ represents the Pauli repulsion 
between the electronic cloud of the molecules and the 
surface electrons, whereas the second term expresses the 
Van der Waals attraction. The values of the parameters 
U , a, Cydm and z v d w used here are the ones supplied in 
Ref. Hi 

The above model is, clearly, highly simplified. By far 
the most important simplification consists of the neglect 
of substrate corrugation, whose role is significant for at- 
tractive substrates such as graphite 20 , but can be ex- 
pected to be less important on substrates such as those 
of Alkali metals, which are relatively weak. 

III. METHODOLOGY 

The Path Integral Monte Carlo method is a numeri- 
cal (Quantum Monte Carlo) technique that allows one to 
obtain accurate estimates of physical averages for quan- 
tum many-body systems at finite temperature. The only 
input of a PIMC calculation is the many-body Hamilto- 
nian Q (which includes the potential energy functions 
U and V). PIMC yields results that are affected only 
by a (small) statistical uncertainty, as well as by an er- 
ror due to the finite size of a many-particle system that 
can be practically simulated on a computer. Comput- 
ing facilities commonly available nowadays comfortably 
afford PIMC computations for systems of several hun- 
dred particles; this size is normally sufficient to obtain 
fairly accurate estimates of energetics and most struc- 
tural properties of interest. PIMC simulations have been 
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extensively adopted to study physisorption of Helium on 
a variety of substrates, including alkali metals 1 ^. 

Because thorough descriptions of the PIMC method 
can be found in the literature 2 !, it will not be reviewed 
here. We only provide some relevant details, namely: 

1. We have experimented in this work with a high- 
temperature approximation for the many-body 
density matrix G(R, R' , r) which is accurate up 
to order t 4 (see, for instance, Ref. l28l) . This 
form is not as effective as the pair-product ap- 
proximation (PPA) , more commonly used in PIMC 
simulations 27 of highly quantal fluids, such as He- 
lium. We have empirically observed convergence of 
the energy estimates with a time step r = (1/320) 
K _1 , almost a three times shorter than in compara- 
ble PIMC studies based on the PPA1&22, The use 
of the high-temperature density matrix considered 
here simplifies (in our view) the calculation, with- 
out rendering it unacceptably inefficient computa- 
tionally. Accurate estimates for quantities other 
than the energy, such as film density profiles, can 
be usually obtained with twice as large a time step. 

2. A multi-level Metropolis scheme based on the stag- 
ing algorithm 2 ^ with free-particle sampling was 
used to sample many-particle paths through config- 
uration space. This technique has been shown to be 
an effective remedy to the inefficiency of single-slice 
sampling, at low temperature. We adjusted the 
length f of the portion of individual single-particle 
paths that the algorithm attempts to update, so 
that the acceptance rate remain above 20%. Typi- 
cally, f w 0.025 K _1 . Rigid displacements of single 
particle paths were performed as well. The average 
displacement was chosen to ensure an acceptance 
rate of approximately 50%. 

3. We have included in our calculation the sampling 
of cyclic permutations of groups of p-B. 2 molecules, 
based on an algorithm very similar to that illus- 
trated in Ref. 27* This procedure is not restricted 
to sampling groups of 2, 3 or 4 molecules, but allows 
in principle cycles of arbitrary lengths (up to TV) to 
be constructed. As mentioned in the introduction, 
however, we found, at the physical conditions ex- 
plored in this work permutations occur exceedingly 
rarely. This could be the result of using less than 
optimal a sampling algorithm, but we have verified 
that our procedure is capable of sampling permu- 
tations, with reasonable efficiency, for bulk 4 He (a 
simulated system of 64 atoms) at T=2 K and sat- 
urated vapor pressure. Thus, we believe the ab- 
sence of permutations in our simulation to reflect 
a physical effect, in agreement with what observed 
by other authors 1 ^* 1 ^. 

Physical quantities of interest, besides the kinetic, po- 
tential and total (e) energy per molecule, are the P-H2 
density profiles as a function of the distance z from the 



substrate, i.e. 

n(z) = i J dx dy p{x, y, z) (3) 

where A is the area of the substrate and p{x, y, z) is the 
three-dimensional (3d) P-H2 density, as well as the the 
angularly averaged, "reduced" pair correlation function 
g{r), with r = \/ x 2 + y 2 and 

s( x ' y) = ~jQ2 J dx ' d y' n ( x + x 'iV + y') n ( x 'i y') ( 4 ) 

where 9 = N/A is the coverage and with n(x, y) = 
J dz p{x,y,z). The correlation function g(r) provides 
a quantitative assessment of the two-dimensional char- 
acter of an adsorbed film; the more two-dimensional an 
adsorbed film, the more closely g(r) mimics the pair cor- 
relation function of a strictly two-dimensional system of 
the same coverage. 

As for the typical size of the physical systems simulated 
in this work, for 6 < 6 = 0.072 A" 2 , PIMC calculations 
were carried out on a system of 36 p-H 2 molecules, ini- 
tially arranged on a triangular lattice at a distance of 3 
A from the substrate. At higher coverage, the initial ar- 
rangement was taken to be a series of successive solid (tri- 
angular) layers, arranged as in a hexagonal close-packed 
lattice, each layer of 2d density o and with an incom- 
plete top layer. Films of up to three layers are studied. 
At the monolayer equilibrium coverage ff e = 0.070 A" 2 , 
we also obtained results for a system of 144 particles. 
No noticeable dependence of the estimates on the size 
of the system can be observed for the physical quanti- 
ties discussed here, with the exception of the energy. In 
all calculations, the height of the simulation box (i.e. the 
distance between the periodically replicated images of the 
substrate) is 30 A, i.e., much greater than the maximum 
film thickness observed, so as to make the use of periodic 
boundary conditions in the z direction (perpendicular to 
the substrate) uninflucnt. 

IV. RESULTS 

PIMC results were obtained in whis work for two sub- 
strates, namely Na and Li. Let us begin with an as- 
sessment of the energetics of the P-H2 films on the two 
substrates. 

Fig. ^ shows computed values of the energy per P-H2 
molecule e(0), as a function of coverage, on a Li and on 
a Na substrate, at a temperature T=2 K. These num- 
bers can be very reliably regarded as ground state esti- 
mates, as values of the energy computed below T w 5 
K are found to be temperature-independent, within the 
statistical uncertainties of the calculation. The dotted 
line in Fig. ^ represents the ground state energy per 
molecule (e& = —85.5 K) of bulk p-H.2, in the fee crystal 
structure, as computed by Silvera and GoldmanSl using 
self-consistent phonon theory^. 

The two curves shown in Fig. ^ share some qualita- 
tive features, e,g., a succession of well-defined, equally 
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FIG. 1: Energy per molecule e (in K) vs. coverage 6 (A -2 ), for 
a p-Eh film on Li (open symbols) and Na (filled symbols) sub- 
strates, computed by PIMC at a temperature T—2 K. Dotted 
line represents the chemical potential of bulk P-H2 at T — 
(from Ref. I30I) . Statistical errors on the PIMC data are of 
the order of the symbol size. 



spaced local minima, corresponding to adlayer comple- 
tion. These minima occur at the same values of 9 on 
both substrates, suggesting that the main player in the 
energy balance is the interaction potential among P-H2 
molecules. 

Let us consider the case of Na first. Although it is (al- 
most monotonically) decreasing, the e(9) curve on Na ap- 
pears to be flattening off, approaching the value eb from 
above, asymptotically. This indicates that no film of any 
thickness will form, i.e., that P-H2 will not wet Na, at low 
T. Clearly, this conclusion can only be made tentatively, 
based on the results obtained up to a value of the cover- 
age 9 = 0.2 A~ 2 . It is possible that a stable (thick) film 
may exist. On performing PIMC calculation on systems 
twice the size of that considered here, one could establish 
more robustly the trend observed here; however, even in 
that case one would need to supplement the information 
provided by PIMC with some other calculation (perhaps 
DF), capable (at least in principle) of making predictions 
in the 9 — > 00 limit. 

Assuming that on Na one indeed has a non-wetting 
situation, one question that arises is that of the temper- 
ature associated to the wetting transition; this issue was 
not addressed in this work. 

On Li, on the other hand, the e{9) curve is monoton- 
ically decreasing for 9 < 9 e , with 9 e w 0.070 A -2 ; it at- 
tains a minimum at 9 e , with e{9 e ) « —100 K, i.e. below 
eh, and displays a slowly increasing trend (while going 
through well-defined local minima, as mentioned above) 
for greater values of 9. It seems physically reasonable to 
expect that e(9) should also flatten off, and approach eb 
from below, in the 9 — > 00 limit, though, again, the finite 
size limitation of this calculation do not allow us to make 
this conclusion firmly. 

The local minima of e(9) correspond to successive solid 
layers (more about this later in the manuscript); be- 
tween these minima, films are thermodynamically unsta- 
ble. This situation is physically quite different from the 
case of a 4 He film on the same substrate, for which con- 



tinuous growth of film thickness as a function of chemical 
potential is observed, with no layering^. 



A. P-H2 monolayer on Li substrate 

In the remainder of this manuscript, we shall concen- 
trate on the substrate on which thermodynamically sta- 
ble P-H2 films are seen to form, i.e., Li. Completion of 
the first layer occurs at a coverage 9 = 0.070 ±0.01. This 
value is consistent with that found on graphite 2 ^, even 
though a graphite substrate is over four times more at- 
tractive, and corrugation plays a significant role. It is 
also only slightly above the equilibrium density^ of pure 
2d p-U 2 . 
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FIG. 2: Typical many-particle configuration for a P-H2 mono- 
layer film adsorbed on a Li substrate (top view) at a temper- 
ature T=2 K. Each black "cloud" represents one of 36 P-H2 
molecules. Periodic boundary conditions are used in both 
directions. The arrangement of the molecules on a triangu- 
lar lattice is clearly seen, as well as the absence of overlap 
among different "clouds", to show that quantum exchanges 
are unimportant in this system, at this temperature. 

This suggests that, as long as the substrate is suffi- 
ciently strong to stabilize a film, the basic physics of 
an adsorbed p-H 2 monolayer film is largely substrate- 
independent, driven primarily by the interaction poten- 
tial between p-H 2 molecules. Fig. [21 shows an instan- 
taneous "snapshot" of a many-particle configuration for 
an adsorbed film at the equilibrium density (top view). 
Each fuzzy "cloud" represents a single molecule; the char- 
acteristic size of a cloud gives a quantitative measure of 
the effect of zero point motion. The arrangement of all 
molecules on a triangular lattice is fairly clear. Also ap- 
parent is the lack of any substantial overlap between dif- 
ferent clouds, to indicate that permutations of particles 
are unimportant in this system. 

More quantitative structural information is furnished 
by the density correlation function g(r), defined in Eq. 
and shown in Fig. [21 The computed g(r) at 9 = 9 e 
is found to be virtually indistinguishable from the pair 
correlation function in strictly two-dimensional P-H2, at 
the same coverage (and at T=0)2i. 

Fig. 01 shows the p-H 2 density profile in the direc- 
tion (z) perpendicular to the substrate (solid line). Also 
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FIG. 3: Reduced pair correlation function g(r), denned 
through Eq. ||IJ, computed by PIMC for a P-H2 monolayer 
film adsorbed on a Li substrate, at coverages 6 = 0.070 A" 2 
and T = 2K. Statistical errors are smaller than symbol sizes. 
The curve is virtually indistinguishable from the pair correla- 
tion function for 2d p-Lh at the same coverage, at T=0. 

shown, for comparison, is the corresponding density pro- 
file for a monolayer 4 He film adsorbed on the same sub- 
strate, at T=0 (from Ref. I33T) . The P-H2 film rests in 
closer proximity to the substrate (by approximately 1.5 
A) with respect to the Helium one, and has about half 
the width. Thus, the physics of a P-H2 monolayer is even 
more 2d than that of one of He. Moreover, the kinetic 
energy per P-H2 molecule at this coverage, in the T — > 
limit is worth approximately 40 K, i.e., ~ 30% more than 
in a purely 2d p-B-2 system at the same coverage^; this 
should be compared to the factor two difference in kinetic 
energy between 2d 4 He and an adsorbed 4 He film on Li 
(from Ref. H3). 

In spite of their lighter mass, P-H2 molecules experi- 
ence a lesser degree of zero-point delocalization in the 
perpendicular direction, because of the much stronger 
attraction to the substrate2&, as well as to the other 
molecules. As more layers are adsorbed on the substrate, 
zero point molecular motion in the perpendicular direc- 
tion is suppressed even further in the underlying layers, 
which become more compressed (see Fig. [5} ■ 



B. Melting of a p-H-2 monolayer 

In order to study the melting of a solid P-H2 mono- 
layer at the determined equilibrium coverage (namely 
8 = 0.070 A -2 ), we computed the energy per molecule of 
the system as a function of temperature (e(T)). Melting 
corresponds to a peak in the specific heat cy = (de/dT), 
The specific heat can show several such "anomalies" , cor- 
responding to different transitions (e.g., solid to solid- 
vapour coexistence); several peaks are indeed observed 
in numerical simulations on graphite, for a commensu- 
rate solid layer with vacancies^. In order to reduce the 
smearing of the peaks, which inevitably occurs when one 
is performing a numerical study on a system of finite 
size, we increased the number of particles from 36 to 144 
in this calculation, corresponding to a simulation cell in 




FIG. 4: Density profile n[z) (A 3 ) in the direction z(A) per- 
pendicular to a Li substrate, for a p-lrb monolayer film of 
coverage 9 e = 0.070 A~ 2 and at T=2 K (solid line). Statis- 
tical errors cannot be seen on the scale of the figure. Dotted 
line shows, for comparison, the corresponding density profile 
for a monolayer 4 He film on the same substrate, at T=0 and 
at the equilibrium coverage 9 e = 0.052 A~ 2 (from Ref. l33f) . 
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FIG. 5: Density profiles n(z) (A -3 ) in the direction z(A) 
perpendicular to a Li substrate, for P-H2 films of varying cov- 
erage, corresponding to one, two and three layers, at T—2 K. 
Statistical errors cannot be seen on the scale of the figure. 

which the area of the substrate A equals 48.74 A x 41.21 
A . Results are shown in Fig. and summarized in Ta- 
ble [I] The energy estimates are not corrected for the 
long-range tail of the interaction potential between P-H2 
molecules, which is truncated at 21 A in the simulation 
(potential truncation is routinely done when studying fi- 
nite systems with periodic boundary conditions). We 
estimate the correction on the total energy to be worth 
no more than 0.03 K. 

The results of Table[l]show that the dependence of the 
total energy per molecule on the temperature reflects es- 
sentially that of the potential energy alone, as the kinetic 
energy depends weakly on temperature, for this adsorbed 
monolayer film. The determination of values for the spe- 
cific heat from the energy estimates can be made by nu- 
merical differentiation, but it is difficult to obtain accu- 
rate results in this way. Moreover, our primary objective 
was to locate the occurrence of the phase transition, and 
the figure clearly suggests a melting temperature T m w 
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FIG. 6: Energy per molecule e(T) for a P-H2 monolayer ad- 
sorbed on a Li substrate. The coverage is 0.070 A" 2 . The 
change of slope of the curve at T ~ 6.5 K is evident. Results 
shown correspond to a system of 144 particles. 
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FIG. 7: Density profiles n(z) (A 3 ) in the direction 2(A) 
perpendicular to a Li substrate, for p-H.2 films of coverage 
8 = 0.070 A -2 , corresponding to one layer, at three differ- 
ent temperatures. The melting temperature is estimated at 
Tm—6.5 K. Statistical errors cannot be seen on the scale of the 
figure. As the temperature is raised, second layer promotion 
of molecules is observed. 



6.5 K, which is where the computed e(T) curve displays 
a sudden, abrupt change of slope. 

These results are similar to those obtained by Wag- 
ner and Ceperley, who studied the melting of a p-H-2 
surfac e 18 ' 19 . The computed melting temperature is also 
quite close to the theoretical 2d result 19 . This suggests 
that this system may provide a rather close experimental 
realization of a 2d system, in many respects even more 
so than an adsorbed Helium film. 

The problem of melting of a 2d solid is interesting and, 
to a considerable degree, still open. In particular, the 
transition from a solid to a liquid phase can be either dis- 
continuous or continuous, the system going, in the latter 
case, through an intermediate, "hexatic" phase featuring 
local six-fold, bond-angle order—. Such a hexatic hase 
can then melt into a disordered fluid as a result of un- 
binding of dislocations 3 ^. The results of this study point 
to a P-H2 monolayer adsorbed on a Li substrate as a can- 
didate to observe such a scenario experimentally. 
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39.40(18) 


-133.79(16) 


-94.39(12) 


7.5 


39.48(10) 


-133.06(17) 


-93.58(12) 


8.0 


39.04(11) 


-132.00(15) 


-92.97(16) 


9.0 


39.16(11) 


-130.45(12) 


-91.30(7) 


10.0 


38.82(8) 


-127.59 (13) 


-88.77(10) 


11.0 


39.43(13) 


-126.76(16) 


-87.33(17) 


12.0 


40.01(14) 


-125.07(13) 


-85.07(13) 


13.0 


40.01(15) 


-122.27(15) 


-82.27(20) 



TABLE I: Kinetic, potential and total energy per P-H2 
molecule (in K) on a Li substrate, computed by PIMC at a 
coverage 6 = 0.070 A -2 at different temperatures. Results are 
for a 144-particle system. Statistical errors (in parentheses) 
are on the last digit (s). These estimates are not corrected for 
the long-range tail of the interaction potential between p-Hb 
molecules, which is cut off at 21 A. We estimate the correction 
on the total energy to be worth no more than 0.03 K 




FIG. 8: Typical many-particle configuration for a p-Eh film 
adsorbed on a Li substrate (top view) at a temperature T—8 
K. Each fuzzy "cloud" represents one of 112 P-H2 molecules in 
the bottom (first) adlayer, the remaining 32 molecules in the 
system having been promoted to the second layer. Periodic 
boundary conditions are used in both directions. The star at 
approximately (29,22) indicates the location of a vacancy. 



In numerical simulations, on the other hand, unam- 
biguous confirmation of dislocation-mediated melting is 
complicated by the finite size of the system that can 
be studied in practice 37 . Fig. shows density profiles 
computed at three different temperatures, below, above, 
and close to the melting temperature. They clearly sug- 
gest that melting occurs concurrently with second layer 
promotion; one can therefore opine that melting occurs 
microscopically with the formation of vacancies in the 
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bottom layer. Visual inspection of instantaneous many- 
particle configurations (an example is shown in Fig. [SJ 
for a temperature T=8 K) at temperatures just above T m 
indeed seems to indicate the presence of vacancies (and 
dislocations) in the bottom layer. 

V. CONCLUSIONS 

Adsorption of P-H2 on Alkali metals substrates was 
studied theoretically by means of Path Integral Monte 
Carlo simulations. In order to assess the expected rele- 
vance of the results presented here to experiments, one 
clearly needs to assess quantitatively the limitations of 
the model utilized, as well those of the calculation. 

The potential used here to describe the interaction be- 
tween two P-H2 molecules, namely the Sil vera- Goldman 
potential, is the one for which most of the previous Monte 
Carlo simulations have been carried out. It is known to 
be reasonably accurate, even though a comprehensive, 
quantitative assessment is, to our knowledge, still lack- 
ing. Its main limitation is that of being merely a two- 
body potential. Although it is designed to incorporate, 
at least in part, effects of three-body interactions, this 
is done only in an effective way. Nonetheless, it seems 
unlikely (to us) that the use of a more refined version of 
the interaction would significantly, qualitatively change 
the results. More important seems the uncertainty of 
the substrate-adatom potential. The potentials used in 
this calculation are the most current, and are a signif- 
icant improvement over the early "3-9" potential, both 
with respect to the functional form as well as the value 
of the most important parameters (e.g., the well depth). 
However, their absolute accuracy has not yet been quan- 
titatively assessed; to that aim, additional experimental 
input may be needed. 

We computed energetics for adsorbed films on both Li 
and Na substrates; results show that thermodynamically 
stable adsorbed P-H2 films will form on Li, whereas a 
Na substrate (and therefore all other Alkali metals sub- 
strates, which are even weaker than Na) is too weakly at- 
tractive for a stable film to form. This result is consistent 
with the experimental observation of a wetting transition 
on Cs and Rb 15 i 16 i 17 . Such a transition is predicted to oc- 



cur on Na and K substrates as well, most likely at a tem- 
perature above the critical temperature, as observed on 
Rb. We also investigated thermodynamic properties of a 
P-H2 monolayer on a Li substrate. This system has an 
equilibrium coverage at low temperature essentially iden- 
tical to that on graphite, approximately 0.070 A~ 2 , and 
is considerably more two-dimensional that a 4 He mono- 
layer adsorbed on the same (Li) substrate, owing to the 
much stronger attraction to the substrate experienced 
by a P-H2 molecule with respect to a 4 He atom. We de- 
termined the melting temperature of the monolayer to 
be approximately 6.5 K, which is very close to theoreti- 
cally computed 2d value. While this value is considerably 
lower than the 3d one (~ 13 K), it is still much higher 
than the temperature at which quantum exchanges are 
expected^ to become significant (~ 1 K). As exchanges, 
and long cycles of particle permutations are known to be 
crucial to the onset of superfluidity^ 7 -, this system does 
not appear to offer a promising avenue for the observa- 
tion of a SF transition in p-H2- Because a Li substrate 
is observed to be relatively strong, whereas a Na one is 
already too weak to allow adsorption, one might wonder 
whether a substrate made of an alloy of Li and Na (or, 
Li and another Alkali metal) might have intermediate 
properties, possibly leading to enhanced zero point mo- 
tion in the perpendicular direction and increased quan- 
tum character. Theoretical investigation of adsorption on 
such a "hybrid" substrate is under way. More generally, 
however, observing novel phases of molecular P-JI2 may 
require achieving a substantial renormalization of the in- 
teraction of p-H 2 molecules, possibly through their inter- 
action with the surface electrons of a metal substrate, or 
of a nanostructure. 
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